00 



< 



13 



WAVE FUNCTIONS AND ENERGY TERMS OF 
THE SCHRODINGER EQUATION 
WITH TWO-CENTER COULOMB PLUS 
HARMONIC OSCILLATOR POTENTIAL 



D.U.MATRASULOV 
Heat Physics Department of the Uzbek Academy of Sciences, 



28 Katartal St., 700135 Tashkent, Uzbekistan 



Abstract 



> 

. Schrodinger equation for two center Coulomb plus harmonic oscillator poten- 

tial is solved by the method of ethalon equation at large intercenter separations. 
Asymptotical expansions for energy term and wave function are obtained in the 



analytical form. 
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a 

The Schrodinger equation with two-center potentials is of considerable 
interest in various problems related with few body systems treating in 
Born-Oppenheimer approximation. It describes the bound states of light 
particle in the field of two heavy particles. Usually such type of systems 
arise in molecular physics. However last years there is a considerable in- 
terest to another systems modelling by two-center Schrodinger equations, 
namely baryons containing two heavy quarks (QQq- baryons)0 and heavy 
flavoured hybryd mesons (QQg-mesons) |5j are becoming subject of exten- 
sive investigation now. The motion of light quark (gluon) in the field of two 
heavy quarks (quark-antiquark pair) can be described in the nonrelativis- 
tic approach by the Schrodinger equation with two-center Coulomb plus 
confinement potential. In this paper we concerned with the Schrodinger 
equation with two-center Coulomb plus harmonic oscillator potential which 



1 



is, to our knowledge, the only Coulomb plus confinement potential allowing 
separation of variables in the Schrodinger equation. Using the method of 
ethalon equation, widely applying for the solution of two-center Coulomb 
Schrodinger equation in molecular physics, we obtain wave functions and 
eigenvalues of the Schrodinger equation with two-center Coulomb plus har- 
monic oscillator potential in the form of asymptotical expansion in terms 
of inverse powers of the intercenter distance. 

So we want to find energy eigenvalues and wave functions of the following 
Schrodinger equation 

[-1a - - - - + u\r\ + rl)]^ = Ei; (1) 
2 n r 2 

In the prolate spheroidal coordinates, which are defined as follows: 
Z = — ^— (K f < oo), rj = —j^— (-K rj < 1) 
this potential can be written in the form 

^ 2 ) = -^f^ + — (2) 



where 



a(f) = 2ZR - ^fe-tt 2 - 1), 60?) = ^ZR - - 1) 

As is well known the Schrodinger equation with potential in the 
form (0) is separable in the prolate spheroidal coordinates. After separa- 
tion we have three ordinary differential equations connected by separation 
constants A and m: 

[|« 2 - 1)| + of + (p 2 - 7 'J 2 )« 2 - 1) - - A]X(f) = 

[|-(„ 2 - i)£ + (p 2 - 7 V)(» 2 - 1) - r^-r, - A]y(u) = o 

or) or) [r) z — 1) 

( ~ +m 2 )Z{(t ) ) 



where 

F 2 



E' = E-^- 1 ' = — 1 a = 2ZR. 

2 ' 4 ' 

Boundary conditions for U and V are 

?/(Ole=i=o, tf(Ol^oo— o, (3) 

V(r]) \ v =±i= 0. (4) 

After substitutions 
these equs. can be reduced to the following canonical form: 

u»(() + £ + ^> - + ^$Mf) = o (5) 

^"(") + [J + 1^2 - *W + (j^Ffo) = (6) 
where a = 2Z/V2E', 7 = w 2 /8£ /2 , 

/i = 2p (7) 

1 Asymptotics of quasi-angular equation. 

We will solve eqs.(|5|) and(g) for large R approximately by the method of 
ethalon equation. This method is sucesfully applied to the solution of 
nonrelativistic two center Coulomb problem |], |2j, || and in the theory of 
difraction of waves. Details of the method of ethalon equation are given in 
@, |3|, |j. Let's start from the angular eq.(g). As an ethalon equation for 
eq.(@) we choose the Whittaker equation 



TT w/ r h 4 h 2 k 1 — m 2 lTir 

W" + [-— + + — rH W 8 

L 4 2 4^ 2 J w 

and seek solution in the form 

V=[z'(rj)]^M k , f (h 2 z), (9) 
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where Mk^(h 2 z) is the solution (regular at zero) of eq.(g). After substitu- 
tion (0) into (0) we get following equation for z: 

4 1{X } /V ( 4 + z 2x(l-x/2) ) + 

^^(i^y-^-i^"^ (10) 

where r = 1 ~! n , x = 1 + 77. Requirement coincidence of trnsition points 

\x=o= 

leads to the following "quantum condition" 

A = 2 ^°) + fi9r 1] 

We will seek the solution of eq. (|H)|) and eigenvalues A in the form of fol- 
lowing asymptotical expansion: 

_ °° zk °° Afe 

fe=0 ^ A:=0 " 

Substitution these expansions into ( |HJ| ) gives us the recurrence system 
of differentil equations for z : 

4 = 2 7 5(x-l) 
= 



1 2kz' {z[) 2 2A 4 2 



2 2 = + 



2^ z 2^ z^x{l-x/2) 2 



and for A 



A = 2*4(0) 
Ai = 2kz' l {0) 

X 2 = 2kz' 2 (0) + M Z M 
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Solving these recurrence equations we have for A 



A (,) = 4kl i + + 0(1) (12) 



and for z 



z = 7 5 X (2 - x) + -L?Zra(l - x) + 0(i) (13) 

hr nr 

From boundary conditions one can obtain for quantum number k |3|, ||j] : 

m + 1 
* = 9 + -2-, 

where g = 0, 1, 2, , .... 

2 Asymptotics of quasi-radial equation 

For equation ([5]) as an ethalon equation we take equation : 

W" + [h 2 s - h 4 y 2 - -^t/V = (14) 



solution of which expressed by the confluent hypergeometric functions [?, ?] 

^ 2 ,s-2c-l 1 

, c H — . 

4 ' 2 



where c = i+vp±l_ 

Boundary condition (g) and properties of functions F |7], §j give rise to 
the following expression for s: 



s = 4n + Vm 2 + 3 + 2. 

Substituting 

U=[y(0]- 1 *W(y(0) 

into eq.([5|) we have 

4r 3-4r^ 1 

= 0. (15) 



h\i 2 -l) 2 4/i 4 7/ 2 2/i 4 ' 
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After substitution 

* 4 

this eq. can be reduced to the form 



where £ = £ — 1 Quantization condition which follows from 4>(x) = is 
written in the form 

A = -2^(0) + I - i[ J + 1] |„ (17) 
Inserting into eq. ( p70|) asymptotical expansions 



°° ^ » A fc 

fc=0 a fe=o " 

and solving obtained equations we have for y: 

y = 2^(t 2 + 2t)2 + -U 7 -^ 2 + 2t)-kn(t + 1) + 



i« 7 -f (* 2 + 2t)-*in?^ + 0(1) (18) 



and for A 



, ( t\ i a 4r — so 7 scry" 3 , , 

A (o = _ 2s7 , _ _ + _ _ _i_ + 0( _) (19) 

3 Asymptotical expansion for energy and wave func- 
tions. 

Asymptotical expansions (\T2j ) and give us expression for energy in the 
form of multipole expansion. In order to obtain this expansion one should 
insert 

E' = E + ^ + ^ + ... 



into (pT2j ) and (|T9j). Equating X 1 and A^ and taking into account (0) we get 
following equations for coeficients E\,E<i,... : 



E x = \[(su - 2ku- 1 )(2E ) 5 i + (4s 2 - 16k 2 - 16t)(2# )*], 
6Z 

E 2 = \e\ + 2su- l E Q + Ei(2E )5Z" 1 (16r 2 + 16A: 2 - 4s 2 ), 



Now we need to find Eq. In order to find this one we note that for 
R — ► oo E' = E Q or 

.2 n>2 

£ = £ + (20) 
On the other hand for large R 

97 oo r E>2 

K(ri, r 2 ) = E( «)^(a>^) + W V + Srtfcostf + — ) + 
R i=o R 4 

n>2 n>2 

(r 2 - 2riteos0 + — )] «w 2 (2r 2 + — ). (21) 
Hence for the energy term with this potential we have 

E = MN + \) + ^f, (22) 



where N = 2n + q + m is the principial quantum number. Comparing ( PD|) 
and (|22|) we have 



E = 2cu(N + ^). 

Thus we have obtained an asymptotical expansion for wave functions 
and energy eigenvalues of the Schrodinger equation with two-center Coulomb 
plus harmonic oscillator potential. Derived formulas have to be useful for 
farther numerical calculations in nonasymptotical region and can be also 
used for estimation QQq baryon energy spectra. 
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